Fate of a Bose-Einstein condensate with attractive interaction 
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We calculate the decay amplitude of a harmonically 
trapped Bose-Einstein condensate with attractive interaction 
via the Feynman path integral. We find that when the num- 
ber of particles is less than a critical number, the condensate 
decays relatively slowly through quantum tunneling. When 
the number exceeds the critical one, a "black hole" opens up 
at the center of the trap, in which density fluctuations be- 
come large due to a negative pressure, and collisional loss will 
drain the particles from the trap. As the black hole is fed by 
tunneling particles, we have a novel system in which quantum 
tunneling serves as a hydrodynamic source. 

1. Introduction 

Bose-Einstein condensation has been observed for the 
alkali atoms 87 Rb 0, 23 Na |], and 7 Li f|, in magnetic 
trapping potentials at temperatures well below 10 _6 K. 
An important difference among these cases is that the 
scattering length between atoms is negative for 7 Li, and 
positive for the others. This indicates a predominantly 
attractive interaction in the former case, and a repulsive 
one for the latter. Thus, we expect the 7 Li condensate 
to be metastable, and will decay with time. The purpose 
of this Letter is to elucidate the mechanism of the decay 
process. 

To be sure, in all three cases the atoms will form 
metallic crystals in free space at low temperatures. This 
means that the short-range interactions are attractive in 
all cases, but the trapping potential tends to keep the 
atoms apart through zero-point motion, and in the case 
of positive scattering length, the atoms do not feel the 
attractive part of the potential. 

We describe the many-boson system with a quantum 
field ip(r, t), which is governed by the Hamiltonian 0j 



H = I d 3 r 



(i) 



where m is the atomic mass, lo is the frequency of the 
external harmonic potential, and u = 47r7i 2 |a|/m charac- 
terizes the interparticle interaction, with negative scat- 
tering length a < 0. The transition amplitude from a 
given initial state ipo(r) to some final state ^i(r) in a 
time interval t\ is given by the Feynman path integral 



T=^ 1 \exp(--Ht 1 ) |Vo> 



(2) 



with boundary conditions ip(r,0) — ipo(r) and ip(r,tx) 
ipi (r). The classical action S[tp] is 



S[lp] = J dt I <l~r 



+ v v v 

at 2m 



— ip ip + -(ip ip) 



(3) 



Physically, tpo(r) and ipi(v) correspond respectively to 
initial and final condensate wave functions, which are 
thermal averages of the quantum field. A Feynman path 
ip(r, t) describes a possible history of the condensate wave 
function. 

We calculate the transition amplitude in the semi- 
classical approximation, in the formal limit Ti — > 0. The 
path integral can be evaluated in the saddle-point ap- 
proximation. The saddle-point path, which corresponds 
to an cxtremum of the classical action, obeys the Gross- 
Pitaevskii (GP) equation with attractive interaction: 



;t .di> 



2m 



-ip — wp*ip 



(4) 



We take as the initial condition ipo(r) = Co exp(— r 2 /2(i 2 ) , 
where Co = \J N/n 3 / 2 d 3 . The normalization is such that 
j d 3 ri/jy = N, where N is the total number of condensate 
bosons. The ground state wave function of the harmonic 
potential has the same form, but with d = do, where do = 

1 /2 

(fi/muj) . For our initial state d is a free parameter. 

2. The tunneling path 

The saddle-point path must be chosen to maximize, 
rather than minimize, the transition probability. Our 
Hamiltonian has no lower bound, and the wave function 
can grow without limit to attain an ever lower energy. If 
the initial wave function is small, however, it must first 
overcome an energy barrier presented by the zero-point 
kinetic term and the external potential. This it can do 



1 



via quantum tunneling, which corresponds to a station- 
ary path in imaginary time We tailor our initial 
discussion to the tunneling process, and put t = —it. 
Self-consistency will tell us whether this is the right thing 
to do. The condensate wave function in imaginary time 
is denoted by <£(r, r) = ^(r,—ir), with complex conju- 
gate <p*(r, t) — ip*(r,ir). The final wave function <fii is 
determined by the time evolution. In other words, we 
follow the path of least action. The transition probabil- 
ity is given by \T\ 2 w A exp (— 2Se/K) 1 where Se is the 
Euclidean action of the saddle-point path: 



Se[4>] = / dr 



d s i 



H Tr~2^ 



(5) 



with T = it. We shall not calculate the prefactor A 
in this paper. Because our Hamiltonian has a trans- 
lation invariance with respect to time, the transition 
probability should be proportional to T (zero mode) as 
T — > oo. Thus the transition rate is proportional to 
exp(— 2H~ 1 Se) to leading order. 

To make all quantities dimcnsionlcss, we measure time 
in units of distance in units of do. Thus, setting 
u> = do = 1, we have the imaginary-time GP equation 



2^ = -V 2 + r 2 0- 5 |^ 



(0) 



where g = 47r|a|/e?o, with initial condition </>o( r ) = 
Cbexp (— ar 2 /2), where a = (d /d) 2 . 

Our problem is quite different from the relativistic case 
[^,0, because our equation is first-order in time and not 
second-order, and hence it describes diffusion rather than 
wave motion. Unlike the relativistic case, the motion 
here does not obey Newtonian mechanics in an inverted 
potential. We therefore do not have either the instanton 
H or the "bounce" @. 

We now make a modified "Thomas-Fermi approxima- 
tion." Noting that initially we have V 2 </>o = (a 2 r 2 — 
3a)(po, we make the replacement V 2 — > (a 2 r 2 — 3a)<fi. 
Assuming that is real, we can now rewrite the equation 
in terms of S 2 : 



ld(f> 2 
2~dr~ 



-V(r)0 2 +g<f> 4 



1 



V(r) = - [(l-a 2 )r 2 +3a] 



The solution is given by 
2 (r,r) 



V(r) 



(7) 



(8) 



g [l + B(r)e 2rf W] 
where B(r) is a monotonically increasing function of r : 

V(r) 1 



B(r) 



9<t>l (r) 



(9) 



For large r we have the asymptotic behavior cf>(r, t) w 

3. Critical number of condensate bosons 

If B(0) > 0, then <j) 2 (r, r) will be positive and finite at 
all times, showing the self-consistency of the imaginary- 
time method. The condition that this be true is 3cn/2 > 
cjCq, or N < N c , with 



JV C = 



8\a\ ' 



The total decay rate is exp (—2Se/K) , where 



Se 



9 Jo 



drr 2 V{r) 



In 



i - m 



-f(r) 



where 



f(r) = 



9<f> 2 a(r) _ N 



exp I 



V(r) N c l + r 2 (l-a 2 )/(3a) 



(10) 



(11) 



(12) 



The critical number of condensate bosons N c has been 
calculated with various methods ]^-[l0|], and all results 
are of the same order of magnitude, and not in disagree- 
ment with experimental indications 111]]. Previous esti- 
mates of the tunneling probability, however, give varied 
physical pictures and mathematical results Jl2|-p^| . Our 
result is closest to that of Ref . Q . 

At each point r in the trap, particles tunnel through 
an energy barrier that depends on r. Consequently, they 
"sweat out" at different rates at different points in space. 
Once over the energy barrier, the saddle-point path veers 
onto the real time axis, and the particles evolve according 
to the PG equation in real time. We shall comment on 
their actual fate later. 

4. The black hole 

When N > N c the initial wave function near the ori- 
gin is large enough to go over the energy barrier. This 
part of the wave function will evolve according to the PG 
equation in real time, and we expect that density fluc- 
tuations will become large, with concomitant collisional 
loss. Thus, a "black hole" opens up near the center of 
the potential, where particles will be drained at a rapid 
rate. In the imaginary-time formalism, the existence of 
the black hole is indicated by the fact that B(r) is neg- 
ative at the origin, and increases to zero at r = b. The 
black hole radius b is determined by the transcendental 
equation 



<t>l(b)=g- 1 V{b) 



(13) 



There is no root unless N > N c . Near the critical value 
we have b 2 w c (N - N c )/N c , where c = 3a/(l + 2a 2 ), 
and for large N, b 2 increases like In AT. 

According to our solution, the condensate wave func- 
tion is independent of r at the black hole radius, de- 
creases with r at any point outside the radius, and in- 
creases with r inside. Inside the black hole, the wave 
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function increases rapidly with decreasing distance, and 
diverges at r = ri (r) , determined by the condition 
1 + £> exp(2rV) = 0. As r increases, ri — > b at an ex- 
ponential speed, so that <f> 2 quickly becomes uniformly 
divergent just inside the rim of the black hole. Further 
into the black hole, </> 2 becomes negative. The last fact 
indicates that the imaginary-time formalism is inconsis- 
tent inside the black hole, and we should go to a real-time 
description. 

5. Negative pressure and instability 

Reverting to real time inside the black hole, we write 
ip(r, t) = y / ne i6 ', and obtain a pair of hydrodynamic-like 
equations One of these is the continuity equation 

dn/dt+V- (nV6>) = 0, and the other reads 09/dt+tt = 0, 
where 



n= l -(yef 



W 777^") ,14! 



Using the modified Thomas-Fermi approximation, we 
replace the second term on the right side by V(r). 
The steady-state solution 80/ dt — gives fi = 0, or 
\v 2 s + V — gn — 0, where v s = X79 is the superfluid 
velocity. This is like Bernoulli's law in hydrodynamics, 
but with negative pressure P — —gn. Solving for the 
superfluid velocity we obtain 



v s = y/2 [gn(r) - V(r)] 



(15) 



The region in which this is real is precisely the black hole. 
Outside this region v s becomes pure imaginary, and we 
are back in the tunneling regime. 

Since the pressure is negative, however, the system 
must be highly unstable with respect to density fluctua- 
tions, which were ignored in our modified Thomas-Fermi 
approximation. The instability must lead to nonuniform 
collapse into pockets of extremely high density, with su- 
perfluid turbulence. We must also take into account the 
fact that inelastic collisions, such as three-body molecular 
recombinations, become increasingly likely with increas- 
ing density, and these collisions will expel particles from 
the trap JlJ,[l7j|. This region is therefore a cauldron of 
density fluctuations and collisional loss, and truly a black 
hole. 

6. Outlook 

The trapped Bose-Einstein condensate with attractive 
interaction offers many exciting possibilities for theoret- 
ical and experimental study. 

(1) The opening and closing of a black-hole region, 
whose size depends on the number of particles, furnishes 
a mechanism for the phenomenon of condensate oscilla- 
tion suggested recently JT^JT^] . A quantitative study will 
necessitate a dynamical treatment. 

(2) In a kinetic theory, we expect that the particles 
"sweating" out throughout the trap, due to tunneling, 



will create a steady flow into the black hole, and drain out 
eventually due to collisions. Here, for the first time, we 
have a system in which macroscopic tunneling serves as 
hydrodynamic source. The kinetic theory in this case will 
have to bridge quantum tunneling and fluid mechanics. 

(3) By plugging the black hole with a laser beam of 
variable size, one can create a "faucet" that controls 
the draining of the condensate. The remaining conden- 
sate, which decays slowly via tunneling, should support 
a larger maximum number. 
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